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ABSTRACT 


The  goal  of  this  thesis  is  to  evaluate  the  use  of  surface  waves  to  detect  buried  objects.  The 
source  used  to  generate  the  surface  waves  was  a  three  element  phased  array  controlled  by 
Lab  VIEW  visual  instruments.  This  research  included  developing  a  source  and  receiver,  evaluating 
attenuation  and  azimuthal  dependence  of  the  surface  wave  propagation,  detection  using  scattering 
within  a  tank  of  sand,  and  using  the  three  element  array  to  beamform  surface  waves.  It  was 
successfully  demonstrated  that  target  localization  using  surface  wave  scattering  and  beamforming 
with  a  phased  array  is  possible. 
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I.  INTRODUCTION 


A.  OBJECTIVE 

The  purpose  of  this  research  was  to  evaluate  the  feasibility  of 
using  a  phased  array  to  generate  surface  waves  to  detect  buried 
land/sea  mines. 

B.  MOTIVATION 

The  detection  of  buried  mines  in  the  ocean  bottom,  surf  zone 
and  on  land  is  a  continuing  civil  and  military  problem.  In  October, 
1994  the  United  Nations  Secretary  General  Boutros  Boutros-Ghali 
commented  that  the  estimated  number  of  mines  still  in  place  from 
conflicts  of  this  century  is  over  100  million  through  within  62 
countries  [Ref  1].  These  mines,  originally  intended  to  disrupt 
commercial  activity,  utility  plant  operations,  port  facilities, 
transportation  routes,  and  troop  movement  take  their  toll  on 
postwar  civilian  populations  throughout  the  world.  The  technology 
to  safely  and  efficiently  remove  land  mines  has  not  been  developed. 
The  current  cost  of  removing  a  single  $25  anti-personnel  mine  is 
estimated  by  the  United  Nations  [Ref  1]  to  ranges  from  $300  - 
$1000.  This  cost  is  projected  to  be  higher  as  popular  plastic  mines 
render  magnetic  mine  clearance  devices  currently  in  use  obsolete. 

The  strategic  advantage  gained  from  the  Iraqi  mine  field  in  the 
Persian  Gulf  can  be  debated,  given  the  results  of  the  war.  However, 
there  is  no  question  that  it  was  effective  in  damaging  two  US. 
warships,  altering  the  assault  objective  to  Faylakah  Island,  and 
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preventing  the  amphibious  assault  on  Kuwait.  At  the  end  of  the 
Persian  Gulf  war  1,300  mines  were  recovered  with  24%  being 
bottom  mines  in  various  states  of  scouring  according  to  the  Center 
for  Security  Strategies  and  Operations  [Ref  2].  These  mines  could 
have  had  a  much  greater  impact  on  the  outcome  of  the  war  and  the 
number  of  personnel  and  material  losses  had  a  different  assault  plan 
been  chosen.  Set  up  as  a  barrier  to  amphibious  assault,  the  Iraqi 
mine  field  was  successfully  avoided  and  removed  once  it  was 
discovered.  If  the  mines  had  been  used  more  effectively  to  hamper 
the  positioning  of  ships  or  if  amphibious  assault  had  been  an  only 
option,  mine  clearance  would  have  been  a  much  bigger  problem  than 
as  it  turned  out. 

As  implied  in  “.  .  .From  the  Sea,”  for  our  current  naval 
doctrine  to  work,  some  method  of  rapid  detection  of  mines  must  be 
developed  to  deal  with  buried  mines. 

C.  BACKGROUND 

For  more  than  25  years  the  US.  Navy  has  utilized  marine 
mammals  to  detect  objects  in  the  water  and  buried  in  the  ocean 
floor.  The  continuing  study  of  these  animals  and  how  their  sensor 
systems  work  could  provide  valuable  insights  for  the  next 
generation  of  sonar  for  mine  detection.  In  a  similar  fashion,  sand 
scorpions  of  the  Mojave  desert,  who  are  able  to  capture  their  prey  by 
detecting  Rayleigh  waves,  might  have  some  of  the  answers  for 
buried  mine  detection,  [Ref  3].  Utilizing  their  eight  independent 
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tarsal  slit  sensilla  as  vibration  receptors,  scorpions  are  able  to 
detect  small  vibrations  in  the  sand  resulting  from  the  Rayleigh 
surface  waves  and  compression  waves  (P  waves)  created  by  the 
movement  of  their  prey.  They  are  then  able  to  accurately  estimate 
time  differences  between  arrivals  of  the  two  types  of  waves  in 
order  to  formulate  target  distance  and  direction  from  the  timing  and 
ordering  of  signals  received  by  their  sensors.  At  distances  of  10  cm 
they  can  derive  target  angle  and  distance  very  accurately  and  tests 
have  shown  they  can  detect  disturbances  out  to  30  cm  [Ref  3]. 

In  elastic  media,  Rayleigh  waves  propagate  with  a  non- 
dispersive  speed  slower  than  the  shear  and  compression  wave 
speeds  of  the  medium.  Because  Rayleigh  waves  are  localized  to 
within  a  few  wavelengths  of  the  surface  of  the  medium  they  could 
be  ideal  for  hunting  mines  that  must  be  placed  near  the  surface  to  be 
effective.  From  1987  to  1992,  BBN  Systems  and  Technology 
undertook  a  feasibility  study  to  design  an  experimental  system  for 
detecting  and  locating  buried  land  mines  using  acoustical  probing 
[Ref  4].  In  their  report,  BBN  employed  source  and  receiver  arrays 
which  achieved  detection  distances  of  15  feet.  The  source  used  in 
that  work  was  a  water  filled  tube  in  which  a  high  voltage  electrode 
created  a  spark.  The  resulting  pressure  wave  was  transmitted 
through  a  rubber  membrane  in  contact  with  the  ground.  It  is 
estimated  that  only  67%  of  the  energy  from  a  point  source  energy  is 
propagated  in  the  form  of  Rayleigh  waves.  The  remaining  energy  is 
disbursed  in  the  form  of  downgoing  vertically  polarized  shear  (SV) 
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waves  (23%)  and  compression  (P)  waves  (10%).  Coupled  with  no 
source  directivity,  BBN's  system  is  inefficient  as  a  method  of 
generating  surface  waves. 

The  focus  of  the  research  presented  in  this  thesis  is  to  build 
an  improved  surface  wave  source.  The  source  will  consist  of  a 
phased  array  of  several  discrete  sources,  operated  at  an  optimal 
frequency  for  the  specific  soil  type. 

D.  STEPS  OF  RESEARCH 

The  stages  of  research  involved: 

1.  Building  an  appropriate  test  tank. 

2.  Developing  sources  and  receivers. 

3.  Develop  data  acquisition  software  using  LabVIEW. 

4.  Measuring  the  performance  of  a  three  element  array  of 
sources. 
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II.  DEVELOPMENT  AND  DESIGN 


The  purpose  of  this  chapter  is  to  1)  review  the  characteristics 
of  Rayleigh  waves,  2)  discuss  the  test  tank,  source  and  receiver  used 
in  this  research,  and  3)  discuss  the  LabVIEW-based  data  acquisition 
system. 

A.  WAVE  PROPAGATION  IN  SOLIDS 

A  point  source  acting  at  the  surface  of  an  elastic  medium 
excites  three  basic  types  of  waves:  compressional  waves,  shear 
waves,  and  surface  waves.  Compressional  waves  are  characterized 
by  particle  motion  along  the  direction  of  propagation,  whereas  in 
shear  waves  the  motion  is  transverse  to  the  direction  of 
propagation.  Both  of  these  types  of  waves  can  exist  within  the  bulk 
of  the  solid  and  spread  spherically  away  from  a  point  source. 

Surface  waves,  as  the  name  implies,  are  confined  to  the  surface  of 
the  medium  [see  Figure  (1.1)].  Their  amplitude  decays  exponentially 
with  distance  from  the  surface,  with  a  decay  constant 
approximately  equal  to  the  wavelength.  As  a  result  of  their 
confinement  to  the  surface,  the  propagation  is  cylindrical  in  nature. 
The  particle  motion  is  a  retrograde  ellipse,  confined  to  the  vertical 
plane.  Figure  (1.2).  The  speed  of  propagation  is  less  than  the  shear 
wave  speed  for  the  medium.  Surface  waves  are  commonly  referred 
to  as  Rayleigh  waves,  although  strickly  speaking,  Rayleigh  wave  is 
the  correct  name  only  for  the  case  of  a  wave  on  the  surface  bounding 
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Figure  (1.1)  Compressional  wave  and  Rayleigh  wave. 


Figure  (1.2).  Rayleigh  wave  propagaton. 


a  semi-infinite  elastic  medium  and  vacuum.  However,  for  practical 
purposes  the  physics  of  Rayleigh  waves  accurately  describes  the 
practical  situation  of  a  wave  on  the  surface  of  a  solid,  many 
surface-wavelengths  thick,  with  air  above. 


B.  TEST  TANK 

In  this  thesis  we  want  to  investigate  the  propagation  of 
surface  waves  at  an  air/wet  sand  interface  in  a  controlled, 
laboratory  environment.  The  first  step  is  to  establish  a  suitable 
environment.  Preliminary  tests  were  performed  in  a  1.1  m  x  3  m, 
fiberglass  lined  wooden  tank,  filled  with  approximately  60  cm  of 
wet  sand.  The  purpose  of  these  preliminary  measurements  were  to 
determine  approximate  surface  wave  speeds,  and  therefore 
wavelengths.  The  results  indicated  wave  speeds  in  the  range  60  - 
90  m/s.  The  corresponding  wavelengths  for  frequencies  from  300  - 
2000  Hz  are  approximately  25-4  cm.  Therefore,  we  need  a  sand 
depth  of  approximately  1  m  to  mimic  a  semi-infinite  medium  at  the 
lowest  frequencies  of  interest.  The  width  of  the  tank  needed  to  be 
as  large  as  possible  to  prevent  interference  with  reflections  from 
its  sides.  These  conditions  were  approximately  met  by  a  twelve 
foot  diameter,  four  foot  deep  redwood  tank.  This  tank  was  filled 
with  approximately  17  tons  of  medium  gain  (.6  mm)  beach  sand  from 
Moss  Landing  in  Monterey  County,  California.  The  sand  was  washed 
to  remove  salt  and  filtered  for  large  objects.  Using  a  garden  hose 
with  spray  attachment  the  surface  was  uniformly  wetted  prior  to 
each  experiment. 

C.  THE  SOURCE 

The  source  design  was  guided  by  the  following  requirements: 

1.  It  should  be  a  vertical  point  source  that  could  generate 
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signals  with  sufficient  amplitude  to  allow  detection  of 
surface  waves  anywhere  in  the  test  tank. 

2.  It  should  provide  a  repeatable  waveform. 

3.  It  should  be  capable  of  high  repetition  rates. 

4.  It  should  have  good  ground  coupling  that  would  last  for  an 
extended  period. 

5.  It  should  provide  clean  signals  over  a  range  of  frequencies, 

from  approximately  300-2000  Hz. 

In  addition,  it  had  to  be  relatively  inexpensive  and  easily 
modified  given  the  time  and  resources  available.  It  was  therefore 
not  the  intent  of  this  research  to  develop  a  working  unit  for  mass 
production.  After  various  types  of  drivers  were  tested  it  was  found 
that  a  four  inch  diameter,  high  cone  excersion,  electrodynamic 
speaker  with  an  aluminum  45  degree  angle  cone  attached  to  the 
center  of  the  diaphragm,  worked  well.  In  addition,  a  25  gram  mass 
was  added  to  the  cone  appendage  to  improve  ground  coupling. 
(Figure(1.3),  factsheet  Appendix  (A). 


Figure  (1.3)  Source  design. 


8 


D.  THE  RECEIVER 


The  design  requirements  for  a  receiver  to  detect  ground 
vibration  for  this  research  were  to  provide: 

1.  Adequate  sensitivity  to  ground  motion  generated  by  the 

source  within  the  distances  of  the  tank. 

2.  Broadband  response. 

3.  Good  ground  coupling  that  could  give  repeatable  results. 

Several  different  designs  were  tried.  The  initial  design 
focused  on  using  a  piezo-electric  disk.  The  disks  were  acquired 
from  small  audio  speakers  available  from  any  commerical 
electronics  shop.  The  disks  had  appendages  glued  to  them  to  provide 
coupling  to  the  ground.  They  also  had  a  mass  attached  with  the  top 
to  improve  coupling  and  reduce  ringing.  Figure  (1.4)  shows  the 
different  appendages  that  were  tried. 


golf  tee 


V 


flat  based 
cylinder 


disk 
resting 
on  surface 


Figure  (1.4)  Receiver  ground  coupling  methods 
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The  golf  tee  allowed  for  the  greatest  sensitivity  however  the 
received  amplitude  varied  by  greater  than  five  percent.  The  inverted 
golf  tee  resting  on  the  surface  with  approximately  20  grams  of  mass 
loading  on  a  ring  element  reduced  variations  to  within  five  percent. 
This  coupled  with  signal  averaging  discussed  later  gave  marginally 
repeatable  results.  However,  the  results  were  not  considered 
satisfactory.  Further  tests  were  done  using  a  commercial  single 
axis  geophone  made  by  Halliburton  Geophysical  Servicves,  model 
number  S4408A.  Though  these  geophones  were  not  as  sensitive  as 
the  various  piezo-electric  sensors,  they  were  found  to  be  more 
reliable  in  producing  repeatable  data. 

E.  SUPPORTING  EQUIPMENT 

The  layout  of  the  test  tank  and  supporting  equipment  is  shown 
in  Figure  (1.5).  The  setup  consisted  of  five  elements: 

1.  The  redwood  tank  filled  with  sand. 

2.  Support  equipment  above  the  tank. 

3.  Lab  instruments  including: 

a.  An  HP  331 4A  Function  Generator  which  was  used  as  a 
source  function  generator  for  attenuation  measurements  and  to 
characterize  cylindrical  spreading.  It  was  also  used  to  provide  a 
reference  signal  when  designing  the  array  function  generator  visual 
instrument  (vi)  created  in  LabVIEW. 

b.  An  COS  6100A  Oscilliscope  to  visually  monitor  analog 
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input  and  output  signals. 

c.  An  HP  3562A  Dynamic  Signal  Analizer  was  utilized  to 
isolate  optimal  propagation  frequencies  for  given  conditions. 

d.  An  Ithaco  1201  Low  Noise  Preamplifier  was  used  to 
filter  and  amplify  to  the  received  signal. 

e.  A  Techron  5507  Power  Supply  amplified  the  source 


signal. 


4.  A  Macintosh  llci  computer  with  a  LabVIEW  interface  board 

provided  data  aquisition. 

5.  Sources  and  receivers 


F.  DATA  AQUISITION  SYSTEM 

It  was  decided  to  perform  data  acquisition  functions  with 
LabVIEW.  The  data  acquisition  functions  are  implemented  through 
software  visual  instruments  or  vis.  The  various  vis  used  in  this 
research  are  explained  later  and  in  the  appendixes.  The  vis  allow 
pulse  generation,  signal  averaging,  filtering,  A/D  conversion,  and  a 
high  data  recording  rate.  In  addition  to  the  software,  two  plug-in 
board  were  needed,  the  NB-MIO-16XL  Data  Aquisition  board  and  the 
NB-AO-6  Analog  Output  board. 


12 


III.  CHARACTERISTICS  OF  A  POINT  SOURCE 


The  first  step  toward  building  a  phased  array  surface  wave 
source  is  to  have  well  characterized  single  sources.  The  purpose  of 
this  chapter  is  to  discuss  the  characterisitics  of  the  waves 
generated  by  a  single  source.  First,  we  examine  the  radial  and 
azimuthal  dependence  the  wave  for  different  frequencies.  Next, 
results  are  presented  of  scattering  measurements  from  a  buried 
target.  Finally,  we  discuss  target  localization  using  multiple 
receiver,  bistatic  scattering. 

A.  RADIAL  DEPENDENCE 

The  purpose  of  measuring  the  radial  dependence  of  the  wave 
propagation  is  threefold:  1)  to  determine  if  the  medium  is  isotropic, 
2)  to  determine  if  the  waves  spread  cylindrically,  and  3)  to 
determine  the  attenuation  coefficient  for  later  use  in  calculating 
beam  patterns  for  the  phased  array. 

1.  Procedure 

The  point  source  was  mounted  on  the  beam  in  the  middle  of  the 
tank  after  the  sand  had  been  wetted  uniformly.  Test  points  were 
laid  out  every  10  cm  along  two  lines,  oriented  at  right  angles.  The 
HP  331 4A  Function  Generator  was  used  to  generate  a  single  cycle 
pulse.  This  pulse  was  amplified  with  the  Techron  5507  Power 
Supply  and  sent  to  the  driver.  The  received  signal  was  sent  to  the 
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the  Ithaco  1201  Low  Noise  Preamplifier  to  be  amplified  and  filtered. 
The  signal  was  then  recorded  using  the  Intensity  Field  vi  as 


shown  in  Figure  (3.3).  Measurements  were  taken  from  the  outer 
radial  distance  moving  toward  the  source  in  order  not  to  disturb  the 
sand  associated  with  successive  measurements.  As  can  be  seen 
from  the  graph,  Figure  (3.1),  the  signal  decreases  monotonically 
with  range.  To  determine  whether  the  spreading  is  cylindrical  and 
the  value  at  the  attenuation  coefficient  (a),  the  measured  amplitude 
is  compared  to  its  assumed  functional  form  [Ref  6]: 
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(1) 


y  =  A! Vr * e 

Dividing  the  amplitude  at  r  by  that  at  a  range  ro  gives: 


y  =  y(r'o)^^lro  I  r^e  (2) 


Rearranging  and  taking  the  natural  logarithm  results  in 


ln(V^  /  ^0*  y  I  yo)  =  -oc(r  -  ro)  (3) 

The  slope  related  to  the  value  of  the  attenuation  coefficient,  has  a 
value  of  2.44  dB/m  which  compares  with  results  for  similar  shear 
waves  derived  by  N.  P.  Chotiros,  [Ref  8],  for  sand,  Figure  (3.3).  By 


Figure  (3.2)  Value  of  attenuation  coefficients. 
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plotting  the  data  in  this  form,  we  can  see  that  data  should  fail  on  a 
straight  line  with  slope  -a.  The  results  for  two  different 
frequencies  are  shown  in  Figure  (3.2). 


Figure  (3.3)  Attenuation  measurements 
by  Chotiros. 

B.  AZIMUTHAL  DEPENDENCE 

Having  investigated  the  radial  dependence  of  the  wave 
amplitude,  the  next  step  is  to  investigate  the  azimuthal  dependence. 

1.  Procedure 

To  insure  relatively  consistent  conditions,  the  tank  was 
watered  and  allowed  to  settle  for  30  minutes  to  improve  moisture 
uniformity.  The  source  was  put  in  place,  and  measurements  were 
taken  as  rapidly  as  possible.  Measurements  were  taken  at  one  meter 
range  with  thirty  degree  seperation  between  locations.  Data  was 
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acquired  using  the  vi  illustrated  in  Figure  (3.4),  and  a  programing 
diagram  in  Appendix  (B).  This  vi  was  designed  to  receive  voltage 
amplitude  triggered  on  the  outgoing  pulse  of  the  HP3314A  function 
generator  and  average  a  specified  number  of  pulses  to  improve 
consistency  in  measurements.  The  results,  for  400  Hz  and  2000  Hz, 
shown  in  Figure  (3.5),  show  that  the  spreading,  was  roughly  uniform 
in  all  directions.  It  is  interesting  to  note  that  these  measurements 
were  taken  utilizing  a  one  cycle  pulse.  Increasing  the  pulse  length 
showed  greater  consistency  in  the  received  signal  amplitude  but 
created  a  larger  problem  with  reflection  off  the  walls  of  the  tank. 
The  results  of  the  measurements,  though  not  ideal,  were  considered 
adequate  to  attempt  beamforming  utilizing  an  array. 

C.  SCATTERING  FROM  A  BURIED  OBJECT 

An  evaluation  of  scattering  from  a  buried  object  was 
conducted  using  a  single  source  and  multiple  receiver  set-up  to 
evaluate  the  feasibility  of  using  an  array.  This  was  accomplished  by 
placing  the  source  at  the  center  of  the  tank  and  three  receivers  at 
approximately  one  meter  ranges  at  various  bearings  out  from  the 
source.  A  vi.  Figure  (3.6),  with  a  diagram  and  explanation  in 
Appendix  (C),  was  designed  to  acquire  the  received  signal  triggered 
on  the  source  pulse  and  average  it  over  a  specified  number  of  pulses. 
The  averaged  received  signal  was  then  saved,  an  a  ten  centimeter 
diameter,  three  centimeters  thick,  copper  cylinder  object  placed  in 
the  tank  covered  by  two  centimeters  of  sand.  The  procedure  was 
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repeated  and  the  second  averaged  signal  was  subtracted  from  the 
first  using  the  vi  in  Figure  (3.7),  with  a  diagram  and  explaination 
included  as  Appendix  (C).  As  can  be  see  in  Figure{3.6),  the  arrival  of 
the  direct  path  and  the  reflected  path  are  revealed  in  the  time 
domain  plot.  By  knowing  the  location  of  the  source  and  the  receiver, 
the  speed  of  the  Rayleigh  wave  can  be  determined  from  the  direct 
path.  The  length  of  the  relected  path  can  then  be  determined.  This 
length  specifies  an  ellipse,  corresponding  to  possible  locations  of 
the  object.  Figure  (3.8).  By  utilizing  multiple  receivers,  an  object 
can  be  located  by  finding  the  intersection  of  multiple  ellipses.  This 
was  accomplished  with  the  data  shown  in  Figures  (3.9,  3.10,  3.11) 
and  plotted  in  Figure  (3.12).  The  location  of  the  source  and  three 
receivers  are  indicated  by  the  open  circles.  The  location  of  the 
buried  object  is  indicated  by  the  x.  It  is  seen  that  the  x  is  near  the 
interections  of  the  ellipses.  This  techique  must  be  developed 
further  but  the  results  are  encouraging. 
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Figure  (3.5)  Azimuthal  dependence  of  waveform 
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Figure  (3.10)  Comparison  of  036.2  degree  path. 
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IV.  PHASED  ARRAY  BEAMFORMING 


Having  investigated  the  single  point  surface  wave  source, 
attention  is  now  turned  to  an  array  of  sources. 


A.  BEAMFORMING 

To  investigate  the  feasibility  of  beamforming  required  some 
changes  to  the  lab  set-up.  A  three-element  array  configuration  was 
selected,  which  required  two  additional  sources  to  be  built.  Because 
the  NB-MIO-16XL  DAQ  card  can  only  support  two  analog  output 
channels,  a  NB-AO-6  analog  output  card  was  acquired  providing 
additional  analog  output  channels.  This  addition  enabled  the  use  of  a 
multichannel  analog  output  control  vi  for  the  array. 


B.  BEAMFORMING  MODEL 

To  predict  the  beam  pattern,  a  model  was  derived 
based  on  the  simple  line  array  [Ref  5].  A  line  array  of  N  simple 
sources  spaced  by  a  distance  d,  having  the  same  amplitude  and  phase 
will  generate  a  velocity  of  the  form 


V{r,  e,t)  = 


j(wt-krt) 


(4) 


where  r/  is  the  distance  from  the  /^^source  to  the  receiver  at 
{r,6).  A  program  written  in  MATLAB,  Appendix  (D),  produced  a  plot 
of  the  beam  pattern  based  on  the  number  of  elements,  their  spacing. 
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amplitude  weighting  and  phase.  The  results  of  this  program  are 
compared  with  the  results  acquired  in  the  sand  tank. 

C.  EXPERIMENT 

The  experiments  were  designed  so  as  to  allow  comparison  with 
the  beamforming  model.  Phasing  and  weighting  for  each  element 
were  identical  in  the  inital  experiment.  To  evaluate  the  ability  to 
beamform,  the  three  elements  of  the  array  were  fixed  on  the  cross 
beam  while  the  receiver  was  free  to  move  azimuthal  at  fixed  range 
from  the  center  element.  The  Beamforming  vi  utilized,  Figure  (4.1), 
diagram  Appendix  (E),  allowed  for  independent  control  of  each 
element  in  frequency,  phase,  amplitude  and  pulse  width  to  a  high 
degree  of  accuracy.  Though  not  utilized  to  its  full  capacity  during 
these  intitial  experiments,  it  is  able  to  support  more  complex  array 
analysis. 

In  the  first  experiment,  the  spacing  was  set  to  one 
wavelength  {X  =  .18  meter)  having  measured  a  corresponding 
surface  wave  speed  (c  =  61 .2  m/s)  at  a  frequency  of  340  Hz.  The 
drive  voltage  and  phase  were  equal  for  all  elements.  The  pulse 
length  was  adjusted  to  fifteen  cycles  to  allow  for  sufficient  overlap 
of  the  signals  from  the  three  sources  to  produce  interference. 

Figure  (4.2)  shows  the  received  waveforms  at  eight  selected 
bearings,  arranged  in  approximately  the  correct  orientation.  The 
center  plot  represents  a  typical  waveform  received  at  bearings 
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3-Element  Radiation  Pattern  for  Theoreticai/Experimental  data 


270 

340  Hz,  15  cycle  pulse 


Figure  (4.3)  Calculated  beam  pattern  for  three  element 
array  overlayed  by  experimental  data. 
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centered  between  the  other  locations.  The  variation  in  amplitude 
from  bearing  to  bearing,  and  the  symmetry  about  the  axis 
demonstrates  the  beamforming.  When  compared  to  the  theoretical 
beam  pattern,  Figure  (4.3)  it  is  seen  that  the  array  is  producing  the 
correct  radiation  pattern.  Figure  (4.4)  shows  the  effect  of 
beamforming  at  a  specific  bearing.  The  upper  figure  is  the  signal 
from  a  single  source,  the  lower  is  the  signal  with  all  three  sources 
driven.  Clearly  beamforming  works. 
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V.  CONCLUSION 


A.  RESULTS 

The  results  of  this  research  are  encouraging.  We  successfully 
demonstrated  target  locatization  using  surface  wave  scattering  and 
beamforming  with  phased  array  technology.  Further  development  of 
the  system  is  recommended. 

The  size  of  the  test  tank  proved  to  be  a  significant  limitation. 
A  next  step  would  be  to  build  a  portable  system  that  could  be  taken 
to  the  beach.  Performing  measurements  at  the  beach  will  allow  the 
use  of  longer  ranges  and  lower  frequencies.  It  will  also  provide 
realistic  background  noise. 

Other  recommendations  are  to  investigate  pulse  shaping  and 
more  advanced  signal  processing,  techniques  to  address  the  issue  of 
probability  of  detection,  etc.  This  project  is  an  ideal  area  to  apply 
recent  advances  involving  time  reversal  signal  processing,  [Ref  6]. 
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APPENDIX  A.  FACTSHEET  ON  SOURCE  ELEMENT 


A  four  inch  woofer  speaker  manufactured  by  Radio  Shack  was  selected 
as  a  source  element.  Keeping  cost  considerations  in  mind,  it  had 
excellent  source  element  characteristics.  The  long  cone  excursion  of 
flexible  rubber  provided  the  necessary  pressure  to  maintain  coupling 
with  the  ground.  The  dust  cap  was  sturdy  enough  to  support  various 
attachments  affixed  without  deforming.  The  backplate  was  a  large 
flat  area  to  allow  for  attachment  of  a  mounting  post,  Figure  (A.1 ). 


WOOFER  SPEAKER 
4'(10cm)  woofer 

FEATURES 

•  Smooth  trequeocy  response 

•  Light  weight  hard  cone  paper  based  on  careful 
study  of  internal  toss  values 

•  Uproiled  rubber  edge  minimizes  distortion  and 
improves  linearity 

•  Long  travel' voice  coil  for  high  power  handling - 


SPECIFICATIONS 

Impedance  :  e  ohms 

Frequency  response  :  50 -7.000Hz 

Free  air  resonance  :  SSHz+SHz 

Ots  .  0.35 

Vas  0.23  eft 

vas  is  the  volume  of  air  that  Is  acoustically 
eqaivalent  to  the  compliance  of  the  cone 
SPL  ;  84+2dB/W  (1m) 

Input  power  nominal)  ;  SW 

Input  power  (maximum)  :  tow 

Magnet  weight  :  228  g  8.0  oz 

Speaker  :  700  g  24.7  oz 


i; 

li 

5 

ii 

;! 

!l 
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HOW  TO  USE 

t.  Installation: 

a  install  the  woofer  unit  from  the  enclosure  front  and  keep  it  flush  with  the  baffle  surface, 
b  Secure  the  speaker  so  that  no  air  leakage  is  permitted. 

2.  Connections  when  using  with  a  tweeter  in  a  two-way  system,  connect  as  shown; 


Note  When  using  the  40-1296  Crossover  choose  the  4000  Hz  crossover  point. 


Figure  (A.1)  Factsheet  on  source  element 
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APPENDIX  B.  LABVIEW  VI  FOR  INTENSITY  FIELD  STUDY 


The  purpose  of  this  vi  is  to  take  field  intensity  measurements 
to  determine  attenuation  and  to  measure  spreading.  The  vi  is 
designed  to  be  triggered  by  the  HP  331 4A  Function  Generator, 
acquire  a  specific  number  of  received  wave  forms  and  average  them. 
The  three  programing  elements  necessary  were  a  "for  loop,"  shift 
registers  and  an  intermediate  vi  call  Al  wave.  The  three 


Graph  or  fiiel 


Figure  (B.1).  Elements  of  Intensity  Field  vi. 


elments  can  be  seen  in  the  simplified  version  above,  Figure  (B.1)  and 
as  elements  of  the  actual  program  in  Figure  (B.2).  The  book  is  a  "for 
loop."  The  arrows  on  each  side  are  the  shift  registers.  The  Al  WAVE 
vi  which  acquires  the  signal  is  also  shown.  In  this  case  the  "for 
loop"  value  is  100.  Therefore,  100  aquisitions  are  made,  summed 
and  then  divided  by  the  "for  loop"  value  to  get  an  average  of  the 
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signal.  The  shift  registers  are  used  to  carry  the  summed  value  to 
the  beginning  of  the  loop  until  the  "for  loop"  value  is  reached.  This 
is  a  timed  acquisition,  meaning  that  a  hardware  clock  is  used  to 
control  the  acquisition  rate.  It  is  also  a  buffered  acquisition, 
meaning  that  the  data  is  stored  in  an  intermediate  memory  buffer 
after  it  is  acquired  from  the  analog  input  channel.  On  the  left  side 
of  the  vi,  in  the  lettered  boxes.  Figure  (B.2),  are  the  front  panel 
inputs  controlling  the  number  of  iterations,  device,  channel  input, 
trigger  type,  and  scan  information. 
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Figure  (B.2)  Field  intensity  vi. 
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APPENDIX  C.  LABVIEW  VI'S  FOR  DATA  ACQUISITION 


Two  vis  were  used  to  acquire  and  analyze  data  to  determine 
object  location  in  the  sand  by  surface  wave  sattering.  The  initial 
aquisition  vi  was  similar  to  the  vi  used  to  take  field  intensity 
measurements  except  it  created  a  binary  file  that  could  be  used  by  a 
comparision  vi.  To  create  this  file  an  addition  was  made  to  the  field 
intensity  measurements  vi  as  shown  below,  Figure  {C.1).  This  can  be 
seen  in  the  vi  on  the  far  right  of  Figure  (C.2). 

The  comparision  vi  was  written  to  read  two  binary  files, 
displaying  the  first  in  the  upper  graph,  both  sets  of  points  in  the 
second  graph,  and  the  difference  of  the  two  in  the  third  graph.  In 
this  way  the  additional  disturbance  in  the  field  could  easily  be 
located  in  the  time  domain  as  shown  in  earlier  Figures 
(3.9,3.10,3.11). 

Once  the  arrival  times  for  the  direct  and  relected  paths  had 
been  determined,  a  plotting  program  written  in  MATLAB  was  created 
to  plot  the  ellipses  with  a  known  distribution.  This  allowed  for 
generating  points  having  the  highest  probability  for  object 
detection. 
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Figure  {C.1).  Elements  of  Data  Acquisition  vi. 


Figure  (C.2)  Data  Acquisition  vi. 
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APPENDIX  D.  MATLAB  PROGRAM  FOR  ELLIPSE 

PLOTTING 

As  shown  in  Figure  (3.8),  it  is  possible  to  create  a  set  of 
points  that  represent  all  the  possible  location  of  the  target  given 
the  arrival  time  of  the  reflected  and  directed  path.  Using  MATLAB,  a 
program  was  developed  to  generate  a  set  of  ellipses  representing 
multiple  receivers.  The  program  relies  on  the  inputs  provided  by 
the  user  shown  seen  in  the  first  few  lines  below: 

%Elipse  Plotting  of  Scattering 

%Create  it  so  the  input  radial  angle  is  the  axis  for  the  elipse 

%Angle  can  be  input  converted  to  x,  y  coordinates 

home 

clear 

clg 

%Receivers 

for  R=1 :3;  %total  number  of  receivers 

bearing  =  input  ('Enter  the  beamg  to  the  Receiver  ’); 
range  =  input  ('Enter  the  range  to  the  Receiver  ’); 

TT  =  input  ('Enter  Travel  Time  of  Reflected  Path  '); 

T  =  input  ('Enter  Travel  Time  of  Direct  Path  '); 

A  =  range*sin((pi/180)*bearing); 

B  =  range*cos((pi/180)*bearing); 

xo  =  A/2;  %center  of  elipse 

yo  =  B/2;  %center  of  elipse 

delT  =  TT-T;  %difference  in  travel  times 

foe  =  range/2;  %center  to  foci  distance 

C  =  range/T;  %wave  speed  based  on  direct  path 

D  =  C*TT;  %distance  traveled  by  reflected  path 

by  =  sqrt  {(DI2)'^2  -  foc'''2);  %minor  axis 

ax  =  D/2;  %major  axis 

t=0:.01:2*pi; 
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a=ax*cos(t);  %ax  is  1/2  length  of  elipse 

b=by*sin(t);  %by  is  l/2width  of  elipse 

X=(a)+xo;  %xo  is  center  in  x  dir 

Y=b;  %yo  is  center  in  y  dir 

Z  =  X+Y*i;  %z=X+Y*i 

Zp(R,:)  =  Z  *  exp(i*(pi*(0.5-(bearing/180)))); 


end 

clc 

dispCPick  Mine  Coordinates  off  the  graph!') 

pause(5) 

figure(l) 

hold  on 

for  R  =  1 :3;  %  Change  for  the  #  of  receivers 

plot(Zp(R,:)) 
end 

plot(0,0,'o') 

plot(76,-6,'o') 

plot(77,10.3,’o') 

plot(52,-38,’o’) 

plot(86,-3i;x') 

hold  off 

text(2,2,’soiirce') 
text(80,-10,  ’#1’) 
text(70,35,  '#2') 
text(50,-50,  ’#3’) 
text(47 ,-22,'object') 
xlabel(’range  (cm)’) 
ylabel('range  (cm)') 
grid 

[K,L]  =  ginput(l) 
clc 

disp(['Mine  is  located  at  X  =  ',num2str(K),'  and  Y  =  ',num2str(L)]) 
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APPENDIX  E.  MATLAB  PROGRAM  FOR  BEAM  FORMING 

ARRAY 


To  accurately  evaluate  the  three  element  array  in  it's  ability 
to  successfully  beamform,  a  MATLAB  program  was  written,  with  the 
help  of  Professor  A.  Atchley,  Naval  Postgraduate  School.  Given 
inputs  of  the  distance  to  the  receiver,  spacing  between  elements, 
wavelength,  and  phase  and  amplitude  for  each  element  it  accurately 
predicts  the  radiation  pattern  of  the  array  in  the  XY  plane.  The 
following  is  the  actual  program  listing: 

%  Plot  program  for  3  element  beam  forming  in  the  far  field 


theta  =  (-pi:.01:pi); 

%angles 

r=l; 

%Distance  to  receiver  from  center 

d  =  .18; 

%Spacing  between  elements 

kd  =  0.9.*2*pi; 

%wave  number  *  d 

k  =  kd/d; 

%wave  number 

%  Source  1  is  at  the  top,  3  is  at  the  bottom. 

fl  =  0*pi/3; 

%phase  control  of  element  1 

f2  =  0*pi/3; 

%phase  control  of  element  2 

f3  =  -0*pi/3; 

%phase  control  of  element  3 

al  =  1; 

%amplitude  of  element  1 

a2  =  l; 

%amplitude  of  element  2 

a3  =  l; 

%amplitude  of  element  3 

%  Range  from  element  to  receiver 

rl  =  sqrt(r*r*cos(theta).*cos(theta)  +  (r*sin(theta)  - 

d).*(r*sin(theta) 

-d)); 

r2  =  r; 
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r3  =  sqrt(r*r*cos(theta).*cos(theta)  +  (r*sin(theta)  + 
d).*(r*sin(theta)  +  d)); 

%  Pressure  created  by  each  element 

pi  =  al  *exp(i*(-k*rl  +  fl))./sqrt(rl); 
p2  =  a2*exp(i*(-k*r2  +  f2))./sqrt(r2); 
p3  =  a3*exp(i*(-k*r3  +  f3))./sqrt(r3); 

pabs  =  abs(pl+p2+p3); 

pnorm  =  pabs./niax(pabs); 

%Experimental  Data 

Rdata=  [380, 110, 110, 110, 130, 130, 150, 150, 380, 380, 280, 
150, 100,... 

150, 150, 50,  80, 180, 360, 280, 100, 100, 120, 100,  80,  80, 

80,  340,... 

280, 140,  no,  130, 160, 130, 140, 150,  380]; 

ttheta=0  :pi/l  8 :2*pi; 
data  =  Rdata/380; 
polar(theta,pnorm); 
hold  on 

xlabel(’375  Hz,  15  cycle  pulse’) 

title(’3-Element  Radiation  Pattern  for  Theoretical/Experimental 

data'); 

grid; 

polar(ttheta,data) 
hold  off 
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APPENDIX  F.  LABVIEW  VI  FOR  BEAMFORMING  ARRAY 


A  LABVIEW  vi  was  written  to  control  the  three  element  array. 
The  basic  design  involved  creating  a  sine  wave  of  one  cycle  with  a 
specific  amplitude  and  frequency  inside  a  "for  loop"  which  generated 
multiple  cycles.  This  was  then  multiplied  by  a  square  wave  that  had 
been  created  in  a  similar  fashion.  The  results  is  a  pulsed  sine  wave. 
These  elements  can  be  seen  in  the  simplified  version,  Figure  (F.1). 


Figure  (F.1).  Elements  of  Three  Element  Function  Generator  vi. 


This  basic  process  is  repeated  three  times,  one  for  each  element 
and  then  combined  into  one  matix  with  each  row  representing  the 
output  of  one  function  generator.  This  matrix  is  then  sent  to  an 
analog  output  write  file  which  repeats  the  data  inside  a  "while  loop" 
until  it  is  told  to  stop.  Figure  (F.2)  is  the  diagram  of  the  working 
program. 
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